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Abstract. Boltzmann's principle S(E, N, V ■ ■ ■ ) = In W(E, N, V ■ ■ ■ ) allows the inter- 
pretation of Statistical Mechanics of a closed system as Pseudo-Riemannian geometry 
in the space of the conserved parameters E, N,V ■ ■ ■ (the conserved mechanical param- 
eters in the language of Ruppeiner [1]) without invoking the thermodynamic limit. The 
topology is controlled by the curvature of S(E, N,V • ••). The most interesting region 
is the region of (wrong) positive maximum curvature, the region of phase-separation. 
This is demonstrated among others for the equilibrium of a typical non-extensive sys- 
tem, a self-gravitating and rotating cloud in a spherical container at various energies 
and angular-momenta. A rich variety of realistic configurations, as single stars, multi- 
star systems, rings and finally gas, are obtained as equilibrium microcanonical phases. 
The global phase diagram, the topology of the curvature, as function of energy and 
angular-momentum is presented. No exotic form of thermodynamics like Tsallis [2,3] 
non-extensive one is necessary. It is further shown that a finite (even mesoscopic) sys- 
tem approaches equilibrium with a change of its entropy AS > (Second Law) even 
when its Poincarre recurrence time is not large. 

1 Introduction 

Why this paper? 

Since more than 100 years does thermo-statistics emphasize the canonical- or 
the grand-canonical ensemble in the thermodynamic limit as the appropriate 
microscopic description of an equilibrized system. Here a homogeneous, practi- 
cally infinite, system is controlled by intensive parameters like the temperature. 
Though in textbooks the microcanonical ensemble is considered as the funda- 
mental ensemble, due to mathematical difficulties it is quickly left in favor of 
the canonical ones. Intensive variables like temperature, pressure and chemical 
potential are used instead of the mechanical defined extensive energy, volume 
and particle number. 

The intensive variables are even canonized to found the axiomatic definition 
of an orthode [4] from where Statistical Mechanics (at least its intensive or 
canonical form) is deduced. Lebowitz [5,6] considers the thermodynamic limit 
and Elliott Lieb [7] extensivity, which also needs the thermodynamic limit to 
ignore surface effects, as the condition sine-qua-non. 
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In the thermodynamic limit limyjy^oo. n/v=p au surface effects may be 
ignored. - Usually - It is clear that this approach cannot do justice to phase sep- 
arations. In fact, the gain in entropy when a system splits into different phases 
by interphase surfaces is the essence of phase transitions of first order. A liquid- 
gas transition is experimentally detected just by the interface between the liquid 
and the gas. Consequently, in the (grand)-canonical approach, phase-transitions 
are indicated by the Yang-Lee singularities [8] where the grand-canonical po- 
tentials are non-analytic in z — or singular. These indicate the break-down 
of the (grand)-canonical formalism. In remarkable contrast, the microcanonical 
ensemble has no problems at phase-separations and the microcanonical density 
of states remains single- valued and multiply diffcrcntiable in all conserved con- 
trol parameters also here, see below. This is certainly the strongest argument in 
favor of the fundamental significance of the micro-ensemble. 

At phase-separation the entropy S(E, N, ■ ■ ■ ) has a positive curvature. Rup- 
peiner's Riemannian geometry of fluctuations [1] must be reformulated there as 
Pseudo-Riemannian. This leads to a negative heat capacity and a violation of 
Clausius' formulation of the Second Law ("heat flows always from the hot to 
the cold system" ) . Phase-separations demand an essential, fundamental, change 
of conventional classical thermo-statistics. Thermodynamics, however, was in- 
vented in the 19. century to explain the working of steam engines. I.e. one of its 
primary issues were just phase-separations. 

What was said applies to large systems with short-range forces. The largest 
systems in nature, self-gravitating astro-physical systems, are subjected to forces 
(gravitation) with a range comparable to the linear extension of the system. 
These systems are naturally inhomogeneous and non-extensive. A description of 
their equilibrium by intensive variables with a homogeneous spatial distribution 
misses these most interesting situations. Ironically, the thermodynamic limit does 
not apply to these really large systems. Tsallis [2] on the other hand addresses 
non-extensive systems explicitely but keeps the description in terms of intensive 
variables which fix the relevant conserved parameters only on average. However, 
non-extensive systems are usually not self-averaging. He believes the equilibrium 
statistics of Hamiltonians systems demands a new definition of entropy. This, is 
not so [9], see section 3.2. 

It is well known that self-gravitating systems collapse to a star in a more 
or less void background at low energies, the "gravothermal catastrophe" [10,11]. 
This is of course quite interesting but outside of any homogeneous canonical 
thermodynamics. There is nothing peculiar with this in microcanonical thermo- 
dynamics. Certainly, angular-momentum can change this essentially. However, 
the microcanonical equilibrium configuration of a self-gravitating system under 
larger angular- momentum has not been investigated yet (exception [12]). I will 
show in section 3.2 how the competition between gravitational collapse and cen- 
trifugal disruption leads in a natural manner to a breaking of rotational symme- 
try and to multi-star configurations with a large variety of different but quite re- 
alistic configurations. In this section the global phase diagram of self-gravitating 
and rotating many-body systems as function of energy and angular-momentum 
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is presented, c.f. [12]. It is for the first time that these various realistic stellar 
configurations are interpreted as global equilibrium configurations. 

The thermodynamic limit is also invoked since Boltzmann to deduce the 
Second Law of Thermodynamics [6,13-15] from microscopic reversible dynamics. 
Then Zermelo's [16] paradox becomes blunted as the Poincarre recurrence time 
t rec is much longer than any physically relevant observation time. This is different 
for a finite, eventually small system. 

By all these reasons a reinvestigation of the microscopic foundation of Statis- 
tical Mechanics starting from Newtonian reversible mechanics of the finite many- 
body system using a minimum of assumptions but avoiding the thermodynamic 
limit is highly needed. The "Geometric Foundation of Thermo-Statistics" pro- 
posed in [17] and further developed here offers a new, deeper, and much simpler 
understanding of the microscopic foundation of Thermodynamics. 

2 The few essentials of Statistical Mechanics 

2.1 Why does conventional statistics like the thermodynamic limit? 

The relative fluctuations of a macroscopic observable A in pure phases of an 
extensive system vanish in the large N limit: 

< A 2 > - < A > 2 1 

We call this behavior self-averaging. Then the probability aspect of statistics 
becomes unimportant. 

Here we want to study also non-extensive situations, therefore, we are not 
allowed to go to the thermodynamic limit. Fluctuations must be taken seriously. 

2.2 Obsolete gospels of conventional thermo-statistics 

Then, most axioms which are mistaken to be fundamental for Statistical Me- 
chanics even turn out to be violated: 

• Phase transitions (do not) exist only in the thermodynamic limit 

• Specific heat is (not) c ~<(SE) 2 > > or dT/dE > 

• Heat does (not) always flow from hot to cold 

• Thermodynamic stability does (not) necessarily imply the concavity of 
S(E,N,---) 

• (No ) extensivity of 5, (no ) scaling with N 

• (No) unique Legendre mapping, T — >E, etc. 

• Rise of entropy is (not) necessarily connected to trend towards uniformiza- 
tion 

• Second Law (not) only in infinite systems, i.e. second law (not) only when 
the recurrence time is much larger than the observation-time t rec 3> t b s 
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• A system at equilibrium is (not) necessarily an orthode in the sense of 
Gallavotti [4]. "Boltzmann's heat theorem", i.e. 

dE(T, P) PdV(T, P) 
T + T 

is (not) necessarily a total differential dS(T, P), because S(T, P) is not always 
a smooth, one valued, function of (T,P), see section 3.1, order parameter. 

Violations of these gospels seem shocking statements: 

Lebowitz [5] and Licb[18,14] believe these make thermo-statistics impossible to 
exist. Nevertheless these violations are valid building stones of statistics. They 
are even necessary for thermo-statistics to apply to the original goal of Ther- 
modynamics, the description of phase-separation, the scenario in which steam 
engines do work. At closer inspection these violations are not so strange. Reca- 
pitulating the main roots of statistical mechanics we will see that it makes a lot 
of sense to formulate it without invoking the above axioms and without using 
the thermodynamic limit. The only essential axiom needed to define equilibrium 
statistics is Boltzmann's principle c.f. section 2.4, eq.(3), once we agree not to 
use the thermodynamic limit. To concentrate on this is a great advantage as this 
principle has a very simple geometrical meaning. 

It is a benefit of the new, extended theory which I offer here, that by reducing 
its axiomatic basis to this single principle it applies also to Hamiltonian non- 
extensive systems and among them to the really large systems as astrophysical 
ones, which are far larger than the thermodynamic "limit" of any homogeneous 
system allows. A whole new world for applications of thermo-statistics opens. 
Of course, then one cannot separate volume from surface effects. This is anyhow 
dubious for non-extensive systems or at phase-separation. 

Here however, I must make it very clear that in any cases where the ther- 
modynamic limit of a homogeneous phase does exist, the geometrical theory 
is in complete agreement with conventional Thermodynamics and conventional 
extensive Statistical Mechanics. 

2.3 Thermodynamics, a probabilistic theory; Control parameters 

Thermodynamics is a macroscopic description of a many-body system within a 
few (M ~ 3) macroscopic control parameters and where (67V — M » M) dof's 
remain uncontrolled. Therefore, Thermodynamics describes all systems with the 
same M simultaneously. All systems under the same macroscopic constraints 
are simultaneously addressed by the theory. Statistical mechanics describes the 
whole 6N — M dimensional manifold, i.e. all points in then iV-body phase-space 
with same energy E , the microcanonical ensemble S 1 . Consequently, it gives 
only probabilistic predictions how the average of the systems in the ensemble 
behave [19]. 

1 We denote manifolds in phase-space by calligraphic letters like £ . 
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A large extensive system in a pure phase is self-averaging. In the thermody- 
namic limit a sharp peak of the probability distribution guarantees the identity 
of the average with the most likely configuration. 

Bur what if the thermodynamic limit does not exist like for a non-extensive 
system? For a small system like a nucleus or an atomic cluster the same mea- 
surement must be performed very often and the average be taken before its 
thermodynamic behavior can be seen. 

2.4 Boltzmann's principle, the microcanonical ensemble 

The key quantity of statistics and thermodynamics is the entropy S. Its most 
fundamental definition is as the logarithm of the area W(E) of the manifold £ 
in the N-body phase-space by Boltzmann's principle [20]: 



Boltzmann's principle is the only axiom necessary for thermo- statistics. With 
it Statistical Mechanics and also Thermodynamics become geometric theories. 
For instance all kinds of phase-transitions are entirely determined by topological 
peculiarities of £(E, N, ■ ■ ■ ) and thus of S(E, N, ■ ■ ■ ) see below. 

Local or global constraints? 

In microcanonical statistics the "extensive" , better conserved, control parame- 
ters as energy, volume, particle number etc. can be determined macroscopically 
sharp. There may well be small, microscopic violations of some microscopic con- 
servation laws due to the non-ideal nature of the container. Therefore, we allow 
small fluctuations in these microscopically conserved quantities. Evidently, it 
does not matter whether the entropy S(E, N, V) has an extremum or not. Its 
local value is significant. It is uniquely defined by eqs. (2,3) as a high-dimensional 
integral. It is thus everywhere multiply differentiable, - certainly the most im- 
portant difference to canonical statistics. This is especially important at phase- 
separations where the curvature of S(E, N, V) is positive c.f. section 3.1. 

Why not canonical? 

In the canonical statistics, also in Tsallis "non-extensive statistics" e.g. the en- 
ergy is fixed in the mean by Lagrange parameters like (3 = 1/T. This works only 
if the microcanonical ensemble is self-averaging. Now for non-extensive situa- 
tions like at phase-separations e.g. the energy E(T) as function of the Lagrange 



W(E,N,V) =e trS(E-H N ) 





S=k-lnW 



(3) 
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parameter like T is multi-valued c.f. section 3.1(ambiguity • • • ). This leads to 
the (in view of standard Thermodynamics) surprising negative heat-capacity c.f. 
section 3.1 (negative heat capacity • • • ), which is well documented even experi- 
mentally c.f. [11, 21-33]. 2 

At points of negative curvature of S(E, TV, • • • ) the canonical ensemble 
is not an orthode 

In chapter 1.5 of his book[4] Gallavotti presented an axiomatic deduction of 
thermo-statistics starting from the definition of an orthode. 

Following Boltzmann's heat theorem [34] he defines an ensemble to be an 
orthode when an infinitesimal change of the energy dE and volume dV makes 

dE_ PdV 

T + T 

an exact differential, at least in the thermodynamic limit. Here T the "temper- 
ature" is the average kinetic energy per particle and P the "pressure" is defined 
as the average momentum transfer per unit time and unit surface area of the 
container. 

Clearly this definition is of little help for a non-extensive system and/or 
when the thermodynamic limit does not exist, where a given T or P docs not 
fix the energy or volume. It fails in situations where the canonical ensemble is 
not equivalent to the microcanonical one, i.e. also at phase-transitions. 

3 Equilibrium Thermodynamics 

3.1 Phase transitions 
Relation to Yang-Lee theory 

In conventional extensive thermodynamics phase transitions are indicated by 
the Yang-Lee zeros of the grand-canonical partition sum Z in the thermody- 
namic limit. In order to generalize the definition of phase transitions also to 
non-extensive systems I start for the moment with the Laplace transform from 
the microcanonical density of states e s to the grand-canonical one (here the 

2 In fact the paper [24] pointed explicitely to the fundamental failure of the canonical 
ensembles near first order phase transitions in general and its non-equivalence to 
the fundamental micro-ensemble which shows a negative heat capacity there. It was 
rejected by Gary S. Grest, the Divisional Associate Editor of statistical mechanics of 
PRL, June 3 1997, with the argument: "I am not convinced that the microcanonical 
ensemble is more physical than the canonical ensemble. After all, phase transitions 
in the experimental world are at constant temperature and not at constant energy." 
If this would be true, a ship would not have been able to sail on the surface of the 
ocean, America would never have been discovered and PRL would not even exist. 



(4) 
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discreteness of N does not matter) : 



Z{T, fl, V) = ir — dN e -\E-»N-TS{E)]/T (5) 

JJo £ o 



V 2 rr°° 

cons.+lin.+quadr. 



■ JJ™ dedne- v ^-» n - Ts( -^' T (6) 



The linear term is explicitely put to by solving 

I - ^ 
T ~ dE 

E — 

T ~ ~ ~dN 

P _ dS_ 
T ~ ~dV 



(7) 
(8) 
(9) 



If s(e,n) is concave (has negative principal curvatures), and there is a single 
solution e s ,n s of cqs. (7) and (8), the stationary point, where integral (6) is 
a double Gaussian integral along the two principal curvatures Vi,v 2 and the 
free-energy density is: 



F(T,»,V) -T\n{Z) + ^N s Tin (y^AQ^/R^) lnV 

V V s s V [ V ' 

The curvature matrix c(e, n) of s(e, n) 

( d 2 s a 2 s \ 

c(e,n) = Jf a at e ( 10 ) 

\ dedn dn' 2 / 

has the eigenvalues :Ai > A 2 ^eigenvectors : Vi,v 2 (11) 

Hessian d(e,n) = \\c(e,n)\\ = AiA 2 (12) 

In general Ai can have either sign. This leads to a new, deeper definition of 
thermodynamic phases: 

Classification by the local topology of curvature 

• A single stable phase is defined by Ai < 0. Here s(e, n) is concave (downwards 
bended) in both directions, fig. (1). There is only a single solution of: 
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Here is a one to one mapping of the (grand)-canonical <->the micro-ensemble. 
The (grand)-canonical has in the thermodynamic limit the same analytical 
properties as the micro-ensemble. It is everywhere smooth, multiply differ- 
entiablc. There arc no discontinuities neither in In Z nor its derivatives. 




Fig. 1. Mono-phase. The order parameter o is defined in subsection "order parameter". 



A transition of first order with phase-separation and surface tension is indi- 
cated by the maximum curvature Ai(e, n) > 0. s(e, n) has a convex intruder 
(upwards bended with Pseudo-Riemannian geometry) in the direction vi of 
the largest curvature (order parameter). There are at least three solutions 
(e s ,n s ): Oi,o 2 ,o 3 see figure (2) of 



1 

T 



[i 



dS_ 

dE 

dS 
dN 



(15) 
(16) 



The whole region {oi, 03} is mapped by eq. (6) into a single point (T tr , v tr ) 
in the canonical ensemble which is consequently non-local in o and bimodal. 
This is the origin of the Yang-Lee singularities. I.e. if the largest curvature 
of s(e, n) is Ai > both ensembles are not equivalent as already pointed 
out by us in 1993[35,36,24,37] see also Barre et al [38]. The possibility of 
positive curvatures (Pseudo Rlemannian geometry) is the main difference 
to the Riemannian geometry proposed by G. Ruppeiner[l,39]. Ruppeiner 
discusses shortly phase transitions of first order and points out that there 
is no divergence of the correlation length £ like at second-order transitions 
and £ remains of the order of the interface thickness, i.e. finite. He does not 
mention (or is not interested in) the positivity of the curvature of s(e,n), 
which encodes important information about the surface tension, see below 
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order parameter 



Fig. 2. Phase-separation, Gibbs double tangent 



and which is also the essential reason why the canonical formalism does not 
apply to non-extensive systems. 
• A continuous ("second order") transition with vanishing surface tension, 
where two neighboring phases become indistinguishable. I.e. where the three 
stationary solutions 01,02, 03 move into one-another. This is at the extremum 
of Ai in the direction of order parameter va=o ■ VAi = 0. These are the 
catastrophes of the Laplace transform eq.(6) E — > T and the critical points 
of the transition. If it is also on the border line (Ai = 0) of the first order 
transition (where Ai < 0), it is the critical end-point of the transition. It is 
an open question whether a line of second-order transition is also the locus 
of the critical end-point of a first-order transition in a hidden parameter [37] . 

Physical origin of positive curvature, the surface tension 

For short-range forces it is linked to the interphase surface tension. This is 
demonstrated for a system of 1000 Na-clusters by figure (3). 
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Fig. 3. MMMC [37] simulation of the entropy s(e) per atom (e in eV per atom) of a 
system of N = 1000 sodium atoms at an external pressure of 1 atrn. At the energy 
e < ei the system is in the pure liquid phase and at e > ez in the pure gas phase, 
of course with fluctuations. The latent heat per atom is qi at = ez — e\. Attention: 
the curve s(e) is artifically sheared by subtracting a linear function 25 + e * 11.5 in 
order to make the convex intruder visible. s(e) is always a steeply monotonic rising 
function. We clearly see the global concave (downwards bending) nature of s(e) and 
its convex intruder. Its depth is the entropy loss due to additional correlations by the 
interfaces. It scales tx N^ 1 ^ 3 . From this one can calculate the surface tension per surface 
atom a sur f /T tr = As sur f * No/N sur f. The double tangent (Gibbs construction) is the 
concave hull of s(e). Its derivative gives the Maxwell line in the caloric curve T(e) at 
T tr , fig- (4). In the thermodynamic limit the intruder would disappear and s(e) would 
approach the double tangent from below. Nevertheless, even there, the probability of 
configurations with phase-separations are suppressed by the (infinitesimal small) factor 

— AT 2 / 3 

e relative to the pure phases and the distribution remains strictly bimodal in the 

canonical ensemble. The region e\ < e < ez of phase separation gets lost. 



Negative heat capacity as signal for a phase transition of first order 

As explained in figure (4) for the example of the q = 10 Potts-model, a positive 
curvature (convex intruder) of S(E) in the energy direction corresponds to a 
characteristic backbending of the caloric curve T(E) or (3(E), and to a negative 
heat-capacity c, the general signal for a phase transition of first order as proposed 
by us more than 15 years ago [41,21,35,25,37,42]: 



dE 




(17) 
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Table 1. Parameters of the liquid-gas transition of small sodium clusters (MMMC- 
calculation [37]) in comparison with the bulk for a rising number No of atoms, N sur f 
is the average number of surface atoms (estimated here as Yl^tluster) °^ a ^ clusters 
with Ni > 2 together. a/T tr = As sur f * No/N sur f corresponds to the surface tension. 
Its bulk value is adjusted to agree with the experimental values of the a s parameter 
which we used in the liquid-drop formula for the binding energies of small clusters, c.f. 
Brcchignac et al. [40], and which are used in this calculation [37] for the individual 
clusters. 





iVo 


200 


1000 


3000 


bulk 




T tr [K] 


940 


990 


1095 


1156 




qiat [eV] 


0.82 


0.91 


0.94 


0.923 


Na 


Sboil 


10.1 


10.7 


9.9 


9.267 




As surf 


0.55 


0.56 


0.44 






N sur f 


39.94 


98.53 


186.6 


oo 




cr/T t r 


2.75 


5.68 


7.07 


7.41 



This was later-on further elaborated by Chomaz and Gulminelly [28,43-46] and 
experimentally verified [30,31]. 

Connecting such a system at e\ + 5E,T\ with another one at 63 — 8E > 
ei + 5E, T 3 and T\ > T 3 then the latter one heats up to T com bi n ed = T3 + AT 
by loosing energy, whereas the former one cools down to T com bined = T\ — AT 
by gaining energy. Here one of the Clausius formulation of the Second Law gets 
invalidated: "Heat flows always from the hotter to the cooler body. " Or with 
other words in the region of negative heat capacity a system acts in equilib- 
rium as a refrigerator. This is well within ordinary classical thermodynamics! 
However, here Ruppeiner's Riemannian geometry [1] must be extended to a 
Pseudo-Riemannian geometry. 
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Fig. 4. a) Specific entropy 

s(e) = / e f3 micro (e)de vs. 
the specific energy e for the 
2-dim. Potts model with 
q — 10 spin orientations per 
lattice point on a 100 * 100 
lattice. In order to visualize 
the anomaly of the entropy 
the linear function a + be 
(a = s(0.2119), b = 1.4185) 
was subtracted. Because 
we use periodic boundary 
conditions one needs two 
cuts to separate the phases 
and the depth of the con- 
vex intruder is twice the 
surface- ent r opy . 

b) Inverse temperature 
/3 m icro{e) = 1/T(e) as di- 
rectly calculated by MMMC 

c) Specific heat 
c(e) = -p 2 /(df3/de). 
The canonical ensemble of 
the bulk jumps over the 
shaded region between the 
vertical lines at ei and e^. 
This is the region of the 
coexistence of two phases 
one with ordered spins, the 
other with disordered spins. 
Here c(e) has two poles 
and in between it becomes 
negative. Canonical ther- 
modynamics is blind to this 
region. Observe that the 
poles are inside e\ < e < e^, 
i.e the canonical specific 
heat (non-dashed region) 
remains finite and positive 
as it should, from [25]. 



The topology of the curvature c(e, n) de-covers the global phase- 
diagram 

Figure (5) shows the example of a micocanonical global phase diagram of the 
Potts (q = 3) lattice-gas as function of the energy e = E/L 2 per lattice-point 
and the relative occupancy n — N/L 2 [47] . L x L is the size of the lattice taken 
to be L = 50, and < TV < I? is the number of occupied sites. 
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The Hamiltonian of the lattice gas is: 



H 



with the occupancy o,; = 



= L~ 2 N = L- 2 J2o t . 

i 

1 , spin particle in site i 



(18) 



, vacancy in site i 



1 

. 8 
. 6 
. 4 
0.2 






1.5 



1 



. 5 







Fig. 5. Global phase diagram or conture plot of the curvature determinant (Hessian), 
eqn. (12), of the 2-dim Potts-3 lattice gas with 50 * 50 lattice points, n is the number of 
particles per lattice point, e is the total energy per lattice point;The line (-2,1) to (0,0) is 
the ground-state of the lattice-gas eq. (18) as function of n. The most right curve is the 
locus of configurations with completely random spin-orientations (maximum entropy). 
The whole physics of the model is between these two boundaries. At the green lines 
the Hessian is det = 0, this is the boundary of the region of phase separation (the 
triangle AP m B) with a negative Hessian. The region of Pseudo-Riemannian geometry; 
At the blue lines is a minimum of det(e, n) in the direction of the largest curvature 
(v\ mam ■ V det = 0), these are lines of second order transition; In the triangle AP m C 
is the pure ordered (solid) phase (det > 0); Above and right of the line CP m B is the 
pure disordered (gas) phase (det > 0); The crossing P m of the boundary lines is a 
multi-critical point. It is also the critical end-point of the region of phase separation 
(det < 0). The red region around the multi-critical point P m corresponds to a flat 
(cylindric) region of det(e,n) ~ and V Ai~ 0, details see [47]; C is the analytically 
known position of the critical point which the ordinary q — 3 Potts model (without 
vacancies) would have in the thermodynamic limit N — > oo, from [47]. 
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Order parameters 

Definition: In the geometric theory, the order parameter o of a phase-transition 
is defined as the length of the trajectory along the direction of maximum cur- 
vature in the global phase-diagram see fig (6). Progressing in that direction one 
experiences the transition from one phase to the other. In the thermodynamic 
limit (when it exists) the region of positive maximum curvature is jumped over 
by the canonical ensemble and the order parameter jumps here also (traditional 
definition of the order parameter). It is important to notice that the order pa- 
rameter is not a simple linear function of the control parameters like E, N c.f. 
fig. (6). 

If there are more control parameters (e,n, ■ ■ •) there might be a situation 
where because of some underlying symmetry the eigenvalue of largest curvature 
Ai is degenerate. In those cases the order parameter is multi-dimensional. All 
these features, convex regions, curvatures of s(e, n, ■ ■ ■ ) are topological properties 
of the entropy surface s(e,n, ■■■) determined by the area of the manifold of 
constant energy etc. within the N-body phase space. 




-2 -1.5 -1 -0.5 

e 

Fig. 6. Various trajectories of maximum curvature Ai, v± =order parameter, for the 
global phase-diagram of fig. (5). Left and below the multi-critical point P m , in the 
region of phase-separation (positive maximum curvature Ai) we see an approximately 
linear behavior. Here /3(e, n) =const. is approximately paralell to z/(e, n) =const. and 
paralell to vi(e,n). 



Ambiguity of intensive parameters and the canonical ensemble 

Fig. (7) explains what happens if one plots the entropy s(e, n) vs. the "intensive" 
quantities (3 = dS/dE and v — dS/dN as one would do for the grand-canonical 
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beta 




Fig. 7. Plot of the entropy s(e, n) as function s m icro{/3, v) of the "intensive" variables 
(y — ds/dn = — f3(i) in the figure labeled as "nu" and /3 = ds/de is called "beta". 
The lines which build the surface are lines for =cons. The positions of the points 
A,D,B,C defined in fig. (5) are only roughly indicated. The convex intruder between 
the lines AP m B and the ground-state line ADB of fig. (5), where s micro (0, v) becomes 
multi-valued as function of v > vp m and > f3p m , here seen from the side, is indicated 
by shadowing. This corresponds to phase-separation at a first order transition. At the 
bottom the projection of the entropy surface onto the {0, v} plane is shown as contour 
plot (lines of equal S(0, v)). The convex part (region of phase-separation) is hidden 
behind the dark "critical" line in the (0, ^)-plane. 



ensemble: As there are several points e^rij with identical 0,v, s m i cro (0,v) is 
a multi- valued function of /?, v . The entropy surface s m i cro (e, n) is folded onto 
itself. In the projection in fig. (7) on the 0, ^-plane, these points show up as a 
black critical line (dense region) . This black line continues over the multi-critical 
point P m towards C indicating the direction towards the critical point of the 
ordinary q = 3 Potts model at n = 1 (zero vacancies). Between P m and C the 
slopes 

vs ne ] (19) 

(20) 

vs ee ] (21) 



ds_ 

00 
or 



ds_ 

dv 



— d ^ Sen 
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are negative large but finite. 

The information given by the projection would be all information which 
can be obtained from the conventional grand-canonical entropy s(T, /i, V), if we 
would have calculated it from the Laplace transform, eq.(6). 

The back folded branches, the convex intruder of s(e, n), the region of phase- 
separation, is jumped over in the Laplace transform eq. (6) and gets consequently 

lost in Z(T,p). Here s(T,p) becomes non- analytical > Yang-Lee singularity. 

This demonstrates the far more detailed insight into phase transitions and crit- 
ical phenomena obtainable by the geometrical interpretation of microcanonical 
thcrmo-statistics [17] but which is not accessible by the canonical treatment, c.f. 
the similar arguments of Gibbs [48]. 

3.2 Rotating self-gravitating systems 
Stars and multi-star clusters 

The most interesting and important non-extensive systems are self-gravitating 
ones. I.e. systems with the Hamiltonian 



with the gravitation <P({ri}) = —Grn 2 J2i<j l r « — r j\~ 1 - r *> Pi an d m denote, 
respectively, the position, momentum and mass of the i-th particle. Because 
of the long range of the gravitation the total potential energy is oc N 2 and 
consequently non-extensive. 

The statistical equilibrium of self-gravitating systems without angular-momen- 
tum was first considered by Thirring [11] who pointed out that microcanonically 
these systems have a negative heat capacity and therefore the microcanonical 
and canonical ensembles are not equivalent. 

Here, we overcome the simple Thirring model and investigate the equilib- 
rium of rotating, self-gravitating systems in a spherical symmetric box under 
various energies E and angular-momenta L. These calculations arc done by 
E.V.Votyakov [12]. 

The following approximations are used: 

a) Mean-field approximation: We approximate any TV-body spatial density pN = 
p(ri, . . . , rjv) by the non-correlated product of single-particle densities [12] and 
work henceforth in the single-particle /x-space. 





JV 



Pn(vi, ■ ■ -,r N ) 



n^) 



(23) 



1=1 



and the gravitation interaction: 



$[Pn] 




(24) 
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The aim is now to find those density profiles p(r) that maximize the (mean- 
field) entropy (k = 1) 

S%(E,L)=]nW%(E,L) (25) 

being the mean- field approximation to the microcanonical "partition sum" , 
i.e. the sum of all uncorrelated many-body states (WJ^ 1 < Wn). After integrating 
over the N -momenta: 



Wj?{E,L) = -£-J E-\L T \^L-<l>[p\ 



P[p]dp(r) (26) 



where P[p] is the probability to observe a density profile p = p(r). 

We use the same trick as Lynden-Bcll [49] and avoid configurations with 
high densities where other physical processes like nuclear reactions become more 
important than gravity. To achieve this, we subdivide the spherical volume V 
into K identical cells labeled by the positions of their centers. The idea is to 
replace the integral over V with a sum over the cells. In order to avoid overlap- 
ping and to cure the short-distance singularity of the Newtonian potential we 
assume that each cell may host up to n particles (1 <C n <C N but Kno > N). 
This condition is essentially equivalent to consider hard spheres instead of point 
particles. The probability P[p] to find a given density profile p(r) is now pro- 
portional to the number of ways in which our N particles can be distributed 
inside the K cells with maximal capacity no and individual occupancies n(r k ) 
or p(r k ) = n(r k )K/V . A simple combinatorial reasoning leads to 

PM = n („;*;,)■ ™ 

colls k,J2 n(r k )=N V V KJ/ 

This looks analogous to a Fcrmi-Dirac statistics, however, here only in coordi- 
nate space. This guarantees the strict non-overlapping condition, p can nowhere 
be larger than p = n K/V. In contrast, in a recent paper Chavanis and Ispola- 
tov [50] use Fermi-Dirac statistics in phase-space. 

We express as a functional of the density profile p(r), such that 

/ p(r)dr = N (28) 
Jv 

and subsequently find the p(r)'s that maximize S^. This leads to the self- 
consistent integral equation: 

log P{X ] =-^U{x) + \p(uy.xY-p (29) 
po-p(x) t> 2 



or, equivalcntly, 



P[X) ~ 1 + e ^(X)4« W xX)H, (30) 



18 D.H.E. Gross 



where uj[p] is the angular velocity (related to the total angular momentum by 
the relation L = I[p]u)[p]), and (3{p\ and U(x) are respectively defined as 

(3= : T (31) 

E-\L T (I\p])-iL-$\p\ 

m = \j U(x)p(x)dx (33) 
b) To solve eq.(30), Votyakov [12] expands p{r,i) into spherical harmonics: 

p(n) =^2bi, m (\n\) Yi m (Gi,<f>i) (34) 

and ignores for simplicity reasons all odd I, m as well also all I > 16. I.e. he allows 
only for parity even, upside - down symmetric configurations (this is later over- 
come). The expansion into spherical harmonics has the advantage that the orig- 
inal non-linear self-consistent three-dimensional integral equation (30) becomes 
now a finite self-consistent set of coupled one-dimensional and two-dimensional 
integral equations. As function of energy and total angular- momentum the mi- 
crocanonical global phase-diagram defined by the topology of the curvature (Hes- 
sian) of S(E, N) shows an astonishing rich picture see fig. (8). 

This example proves the superiority of the geometrical, topological, inter- 
pretation of thermo-statistics. It reproduces realistically many different config- 
urations of even gravitating systems. No canonical description nor Tsallis non- 
extensive statistics had achieved this. 
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Fig. 8. Global phase diagram of a self-gravitating system as function of energy E and 
angular-momentum L in dimension-less units. Systematic calculations were done only 
in the dotted region. There are 3 different mono-phases (regions of negative maximum 
curvature Ai < 0): A single star plus some gas at small E and L; at high energy and 
moderate L, a gas phase; and at large L configurations of broken spherical symmetries 
with double-star systems. The mixed phase is the region with positive maximum cur- 
vature Ai > where at least two different phases (single star and gas (low L), single 
star and double star systems, eventually rings, at intermediate L are in competition, 
finally double-stars and gas compete with one-another at large L and E. 



4 Approach to equilibrium, Second Law 
4.1 Zermelo's paradox 

When Zermelo [16] argued against Boltzmann, that following Poncarre any 
many-body system must return after the Poincarre recurrence time t rec and 
consequently its entropy cannot grow indefinitely, Boltzmann [34] answered that 
for any macroscopic system t rec is of several orders of magnitude larger than 
the age of the universe, c.f. Gallavotti [4]. Still today, this is the answer given 
when the Second Law is to be proven microscopically, c.f. [51]. Then, Zermelo's 
paradox becomes blunted. 

Here, I argue, even a small system approaches equilibrium with a rise of its 
entropy AS > under quite general conditions. Thus, Zermelo's objection must 
be considered much more seriously. 

However, care must be taken, Boltzmann's definition of entropy eq.(3) is 
only for systems at equilibrium. To be precise: in the following I will consider 
the equilibrium manifold £(E, V a ) at t < to. At to the macroscopic constraint V a 
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is quickly removed e.g. a piston pulled quickly out to V a + V b , and the ensem- 
ble is followed in phase-space how it approaches the new equilibrium manifold 
E(E,V a + V b )eeeSg.(9). 

4.2 The solution 

Entropy does not refer to a single point in iV-body phase space but to the whole 
ensemble £ of points. It is the ln(W) of the geometrical size W of the ensemble. 
Every trajectory starting at different points in the initial manifold A4(t = to) = 
£(E,V a ) <G £(E,V a + 14) spreads in a non-crossing manner over the available 
phase-space £(E, V a + Vb) but returns after t rec . Different points of the manifold, 
or trajectories, have different t rec which are normally incommensurable. I.e. the 
manifold A4(t) spreads irreversibly over £(E, V a + V b ). 



When the system is dynamically mixing then the manifold M.(i) will "fill" the 
new ensemble £(E, V a + V b ). Though at finite times the manifold remains com- 
pact due to Liouville and keeps the volume W(E, V a ), but as already argued by 
Gibbs [52,53] Ai(t) will be filamented like ink in water and will approach any 
point of £(E, V a + Vb) arbitrarily close . Then, lim ^oo M(t) becomes dense in the 
new, larger £ (E, V a + Vb). The closure M.(t = oo) becomes equal to £(E, V a + Vb). 
I.e. the entropy S(t = oo) > S(to)- This is the Second Law for a finite system. 

Macroscopic resolution, fractal distributions and closure [19] 

We calculate the closure of the ensemble by box counting [54] . Here the phase- 
space is divided in Ng equal boxes of volume S 6N . The number of boxes which 
overlap with M (t) is Ng and the box-counting volume is then: 



The box-counting method is illustrated in fig. (9). The important aspect of the 
box-counting volume of a manifold is that it is equal to the volume of its closure. 

At finite times A4(t) is compact. Its volume W(t) equals that of its closure = 
W(to). However, calculated with finite resolution 5 > 0, Wg(t) becomes > W(t) 
for t larger than some tg, where 



Mixing 



Q d (S) = N s S d , 

Q d = lim Q d (5). 



here with d = 6N — 1 



(35) 
(36) 




5(E - H N {q(t),p(t); [q(t ), P (t ) € £(V a )}})(37) 



A natural finite resolution would be 
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Of course the actual problem will often allow a much coarser resolution because 
of the insensitivity of the usual macroscopic observables. Then the equilibration 
time t$ will also be much shorter. 

Thus the new definition of Boltzmann's principle eqs.(2) is: 

S = ln(W 4 ), (39) 
or mathematically correct, though unphysical, at infinite times: 

S = lim lim ln(W s (t)) = N\n[(V a + V b )/V a ] + S(t ). (40) 



V a V b V a + V b 













































































t < t t> to 

Fig. 9. The compact set £ (to), left side, develops into an increasingly folded but non- 
crossing "spaghetti" -like distribution £(t, to) in the phase-space with rising time t after 
opening the volume Vb- The right figure shows only the early form of the distribution. 
At much later times it will become more and more fractal and finally dense in the 
new larger phase space. The grid illustrates the boxes of the box-counting method. 
All boxes which overlap with M(t,to) contribute to the box-counting volume and are 
shaded gray. Their number is Ng 



5 Conclusion 

The geometric interpretation of classical equilibrium Statistical Mechanics by 
Boltzmann's principle offers an extension also to the equilibrium of non-extensive 
systems. 

Because microcanonical Thermodynamics as a macroscopic theory controls 
the system by a few, usually conserved, macroscopic parameters like energy, par- 
ticle number, etc. it is an intrinsically probabilistic theory, ft describes all systems 
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with the same control-parameters simultaneously. If we take this seriously and 
avoid the so called thermodynamic limit Qim.v-nx>,N/v=p)i the theory can be 
applied to the really large, usually inhomogeneous, self-gravitating systems. In 
chapter (3.2) it is shown how this new approach enables to view many realistic 
astro-physical configurations as equilibrium configurations under the control of 
total energy and angular-momentum, c.f.[12]. 

Within the new, extended, formalism several principles of traditional Sta- 
tistical Mechanics turn out to be violated and obsolete. E.g. we saw that at 
phase-separation heat (energy) can flow from cold to hot. Or phase-transitions 
can be classified unambiguously in astonishingly small systems. These are by 
no way exotic and wrong conclusions. On the contrary, many experiments have 
shown their validity. I believe this approach gives a much deeper insight into 
the way how many-body systems organize themselves than any canonical statis- 
tics is able to. The thermodynamic limit clouds the most interesting region of 
Thermodynamics, the region of inhomogencous phase-separation. 

Because of the only one underlying axiom, Boltzmann's principle eq.(3), the 
geometric interpretation keeps statistics most close to Mechanics and, therefore, 
is most transparent. The Second Law (AS > 0) is shown to be valid in closed, 
small systems under quite general dynamical conditions. 
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